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A.  Description  and  Motivation  of  the  Scientific  Research  Goals 

To  advance  the  fiindamental  understanding  of  the  electronic  bonding,  atomic  geometry, 
and  vibrational  properties  of  ferroelectrics,  we  undertook  to  develop  a  first  principles  linear 
response  approach  based  on  the  linearized  augmented  plane- wave  (LAPW)  method.  The 
following  goals  during  this  grant  period  were  all  reached.  The  motivation  for  these  goals  and  a 
more  detailed  description  follows. 

1)  Develop  a  first  principles  LAPW  linear  response  method  within 

the  framework  of  density  functional  theory. 

2)  Implement  the  method  into  a  computer  code  that  is  accurate,  efficient, 

and  conveniently  interfaces  with  our  existing  LAPW  bandstructure 

computer  code. 

3)  Benchmark  the  method  on  well-studied  systems  to  establish  the  accuracy  of  the 

method. 

4)  Apply  the  method  to  new  systems. 

Why  use  linear  response  theory? 

All  previous  first  principles  density  functional  studies  of  ferroelectrics  by  our  group  and  by 
others  were  based  on  conventional  supercell  treatments.  These  could  examine  only  limited  aspects 
of  the  phonon  spectrum  and  the  structural  and  vibrational  instabilities  related  to  ferroelectric 
phase  transitions.  In  particular  they  were  limited  to  studying  coherent  and  periodic  structural 
rearrangements  that  could  be  modelled  by  small  supercells.  There  had  been  no  detailed 
examinations  of  the  structural  instabilities  for  general  wavevectors,  although  this  may  be  needed 
to  describe  the  character  of  the  ferroelectric  phase  transition.  The  LAPW  linear  response  method 
that  we  have  developed  during  this  grant  period  is  very  powerful,  because  it  permits  calculations 
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of  the  electronic  response  to  periodic  perturbations  of  arbitrary  wavelength,  while  keeping  the 
computational  burden  similar  to  that  of  a  self-consistent  calculation  for  the  unperturbed  system. 
Thus,  for  example,  phonon  frequencies  may  be  efficiently  mapped  at  all  wavevectors  throughout 
the  Brillouin  zone.  Typically  in  supercell  calculations,  only  selected  phonons  at  the  zone  center 
(infinite  wavelength)  and  a  few  zone  boundary  phonons  could  be  examined.  In  addition.  Bom 
effective  charges  and  the  high-frequency  dielectric  constants,  6„ ,  which  are  difficult  to  obtain  by 
other  methods  are  also  easily  obtained. 

Why  base  the  linear  response  approach  on  the  LAPW  method? 

The  density  functional  linear  response  method  is  relatively  new,  and  most  applications 
have  used  the  pseudopotential  method  with  a  plane-wave  basis.  Ferroelectric  materials  are  very 
challenging  to  a  plane- wave  basis,  however,  because  of  the  presence  of  atomic  species  such  as 
transition  metal  atoms  (e.g.,  niobium  and  titanium)  and  first  row  atoms  (e.g.,  oxygen),  which 
have  localized  bonding  orbitals.  Numerical  convergence  with  the  size  of  the  plane-wave  basis  is  a 
difficulty  for  these  materials.  The  primary  goal  of  this  grant  was  to  implement  the  linear  response 
theory  using  a  more  practical  set  of  basis  functions  without  sacrificing  any  accuracy.  The  LAPW 
method,^  which  has  been  applied  with  great  success  to  a  variety  of  physically  interesting  systems, 
was  chosen  for  this  purpose.  The  LAPW  method  uses  basis  fiinctions  that  are  a  combination  of 
plane-waves  and  localized  functions,  and  it  has  no  difficulty  treating  such  localized  orbitals. 


4 

B.  Significant  Results 

Development  of  the  formalism  and  testing  of  the  computer  code 

During  this  grant  period,  the  formalism^  for  implementing  the  linear  response  method  in 
the  LAPW  method  was  developed.  Developing  the  formalism  was  a  major  task,  since  unlike 
plane-wave  based  methods,  the  LAPW  method  employs  position  dependent  basis  functions.  In 
analogy  with  the  occurrence  of  Pulay-type  corrections  in  the  calculation  of  forces  using  the 
conventional  LAPW  method,^  this  positional  dependence  introduces  additional  quantities  that 
must  be  calculated  in  the  LAPW  linear  response  method,^  significantly  modifying  the  plane- wave 
based  linear  response  formalism. Once  the  details  of  this  new  formalism  were  worked  out  on 
paper,  a  computer  code  was  developed.  It  has  been  thoroughly  tested  by  computing  phonon 
modes  for  wave  vectors  q  throughout  the  Brillouin  zone  in  Si  and  GaAs.^  Detailed  comparisons 
with  conventional  LAPW  supercell  calculations  and  with  experimental  vibrational  frequencies 
confirm  that  accurate  results  are  obtained  and  that  the  method  can  now  be  used  with  confidence. 
Effective  charges  and  the  high-frequency  dielectric  constants,  6„ ,  which  are  difficult  to  obtain  by 
any  other  method,  were  also  calculated  to  be  in  excellent  agreement  with  previous  calculations 
and  experiment  for  these  materials.^ 

Application  to  the  search  for  multi-well  Bom  Oppenheimer  potentials  in  CuCl 

CuCl  is  known  to  exhibit  large  anharmonic  effects  and  possibly  a  complicated  multi-well 
Bom  Oppenheimer  surface  reminiscent  of  the  instabilities  in  perovskite  ferroelectrics,  and  this  was 
the  principle  motivating  factor  for  studying  this  material.  Unlike  conventional  semiconductors 
such  as  Si  and  GaAs,  which  can  also  be  treated  by  the  plane-wave  based  linear  response  approach. 
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a  complete  treatment  of  CuCl,^  has  not  been  previously  published  because  the  localized  3d- 
orbitals  of  the  Cu  make  this  material  difficult  to  treat  with  plane-wave  based  methods.  We 
calculated  its  phonon  dispersion  and  found  it  to  be  in  excellent  agreement  with  experiment, 
despite  the  other  observed  anomalous  properties.  There  was  no  indication  of  multi-well  potentials. 
The  calculated  Bom  effective  charges  were  within  10%  of  experiment  and  the  high  frequency 
dielectric  constant  e.  was  about  30%  larger  than  the  experimental  value  (similar  overestimates  of 
by  density  functional  calculations  are  found  in  simpler  materials  too).  The  overall  level  of 
agreement  with  experiment  is  comparable  to  that  attained  by  other  density  fimctional  calculations 
in  simpler  semiconductors.  Despite  the  anomalous  properties  of  CuCl  and  the  presence  of  large 
anharmonic  interactions,  linear  response  calculations  provide  accurate  harmonic  phonon 
dispersion  curves.  A  similar  conclusion  was  reached  in  the  high-T^  superconductors  LaCuO  and 
YBCO  where  calculated  frequencies  (using  a  conventional  supercell  approach)  are  generally  in 
good  agreement  with  experiment  despite  perhaps  even  greater  anharmonicity.’  This  is  very 
encouraging  for  applying  the  LAPW  linear  response  method  to  perovskite  ferroelectrics. 

C.  Impact  of  this  Work  on  Future  Research 

Since  the  grant  period  ended,  the  LAPW  linear  response  method  has  subsequently  been 
applied  to  study  the  ferroelectric  phase  transitions  in  perovskite  ferroelectrics.  A  recently 
completed  project*  examined  KlSlbOj,  one  of  the  most  extensively  studied  systems  of  the 
perovskite  family  of  ferroelectric  materials.  The  character  of  its  ferroelectric  phase  transitions  (and 
those  of  related  systems  like  BaTiO,)  has  continued  to  be  a  subject  of  controversy.  In  recent 
years,  the  weight  of  evidence  had  shifted  from  a  soft-phonon  description  to  support  of  order- 
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disorder  models,  although  recent  dif&se  x-ray  scattering  experiments  seem  to  support  the  soft- 
phonon  picture.  Previous  density  functional  total  energy  calculations  found  q  =  0  instabilities 
providing  qualitative  support  for  the  order-disorder  scenario.  We  have  discovered  that  the 
instabilities  ofKNb03  (and  most  likely  those  of  related  perovskite  ferroelectrics)  are  considerably 
more  complex  than  previously  suspected.  We  calculated  the  lattice  dynamics,  mapping  out  the 
instabilities  at  all  wavevectors  in  the  Brillouin  zone.  The  principle  results  found  in  the  ideal 
KNbOs  perovskite  structure  are  (i)  the  system  is  stable  against  displacements  of  a  single  Nb  atom, 
and  (ii)  unstable  TO  modes  exist  for  q  perpendicular  to  the  <100>  directions  in  the  BZ.  The 
extent  of  these  modes  in  the  BZ  implies  instabilities  against  formation  of  chains  of  Nb  atoms 
(mutually  uncorrelated)  along  <100>  crystallographic  directions.  The  existence  of  such  chain 
instabilities  had  been  suggested  based  on  early  x-ray  measurements,  but  this  is  the  first  time  that 
first  principles  calculations  have  provided  evidence  for  this  unusual  feature. 
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